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We investigate the strong-coupling regime of the stationary Kardar-Parisi-Zhang equation for 
interfaces growing on a substrate of dimension d = 1, 2, and 3 using a nonperturbative renormaliza- 
tion group (NPRG) approach. We compute critical exponents, correlation and response functions, 
extract the related scaling functions and calculate universal amplitude ratios. We work with a 
simplified implementation of the second-order (in the response field) approximation proposed in a 
previous work [PRE 84, 061128 (2011) and Erratum 86, 019904 (2012)], which greatly simplifies the 
frequency sector of the NPRG flow equations, while keeping a nontrivial frequency dependence for 
the 2-point functions. The one-dimensional scaling function obtained within this approach compares 
very accurately with the scaling function obtained from the full second-order NPRG equations and 
with the exact scaling function. Furthermore, the approach is easily applicable to higher dimensions 
and we provide scaling functions and amplitude ratios in d = 2 and d = 3. We argue that our ansatz 
is reliable up to d ~ 3.5. 
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I. INTRODUCTION 



In their seminal paper Kardar, Parisi and Zhang 
(KPZ) proposed the nonlinear Langevin equation 



dh(t,: 
dt 



= vV 2 h(t, 



+ V(t,x) (1) 



as a continuum model for the dynamics of a growing in- 
terface. In this equation, h(t, x) is a single valued height 
profile depending on the cZ-dimcnsional substrate coor- 
dinate x and on time t. The term r](t,x) represents an 
uncorrelated white noise with zero mean (j)(t,x)) = 
and strength D, 

x)7?(t',f )) = 2DS d (x- x")S(t~t'), (2) 

which accounts for the randomness in the growing mecha- 
nism. The term v~S/ 2 h describes the surface tension which 
tends to smoothen the interface. The nonlinear term, 
which represents an excess growth along the local normal 
to the surface, is the essential ingredient to capture the 
nontrivial behavior of the growing interface. The KPZ 
equation generates generic scaling (without fine-tuning 
any parameter). For instance, the 2-point, height-height 
correlation function in the co-moving frame 



C{t,x) = (h(t,x)h(0,0)) c 



(3) 



endows at large time t and distance x = \x\, the dynam- 
ical scaling form 



C(t,x) = x 2x F(t/x z ), 



(4) 



where \ an( i z are the universal roughness and dynam- 
ical critical exponents and F(y) is a universal scaling 
function. The KPZ equation ([1]) encompasses a nonequi- 
librium phase transition in dimension d > 2 between a 



smooth (x = 0) and a rough (x > 0) phase, separated 
by a critical value A c of the non-linearity. In dimension 
d < 2, the interface always roughens. 

In fact, the KPZ equation arises in connection with 
a large class of nonequilibrium and disordered systems 
and has thus emerged as one of the fundamental theo- 
retical models for nonequilibrium phase transitions and 
scaling phenomena 0, [3| ■ Despite its apparent simplic- 
ity, the KPZ equation has resisted most of the theoretical 
attempts to provide a complete description of its strong 
coupling phase in generic dimensions. Very recently, an 
impressive breakthrough has been achieved to describe 
one-dimensional interfaces, bringing out both exact an- 
alytical solutions and high-precision experimental 
measurements [To|]. From the theoretical point of view, 
the one-dimensional problem is special and much sim- 
pler due to the existence of an additional generalized 
fluctuation-dissipation theorem (FDT) in this dimension. 
Thereby, the exponents are known exactly as z = 3/2 
and x = 1/2- However, in dimension d > 1, the KPZ 
equation remains unsolved. This deceptive situation has 
persisted due to a lack of controlled analytical tools to 
describe the strong-coupling regime. Notably, standard 
perturbative expansions and renormalization group (RG) 
methods (II. [Tll - ll3| were proved to fail at all orders to find 
a strong coupling fixed point fldl ] . Some nonperturbative 
approaches have been devised, such as the mode-coupling 
(MC) approximation [l5rll8j . the self-consistent expan- 
sion (SCE) [lH , or the weak noise scheme [2(1 [2l| . Beside 
critical exponents, they provided some results concerning 
the one-dimensional scaling functions and various (con- 
flicting) predictions regarding the existence of an upper 
critical dimension, above which x = 0, that is the inter- 
face would remain smooth for all values of the nonlinear- 
ity. 

In the present work, we concentrate on the universal 
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scaling functions and amplitude ratios related to the KPZ 
equation. The existence in generic dimensions of a scal- 
ing form for the stationary correlation function of the 
KPZ growth and of universal amplitude ratios was early 
observed [12-113] • This observation was corroborated in 
d = 1 with reasonable agreement by many different ap- 
proaches, including Monte-Carlo simulations of discrete 
models, direct numerical integration of the MC equa- 
tions [l|| [H, H|[ and a soliton approximation [20| ■ In 
subsequent studies, a stretched exponential decay of the 
tail of the scaling function was predicted independently 
by the SCE [27| and the MC [H, HU approximations, 
and numerically verified in one dimension within both 
approaches 0, [3C|. Later, the one-dimensional scaling 
function was calculated exactly in the long-time limit [4j 
(and for all the scaling regime very recently 0]). The 
exact stationary scaling function in d = 1 shows an ex- 
ponential tail, but on a different scale than that brought 
out by MCT and SCE approximations [U, [H, HI- On 
the other hand, in higher dimensions, apart from the pre- 
diction of the asymptotics of the tail from SCE and MCT, 
theoretical calculations of scaling functions or amplitude 
ratios are scarce, if not incxistent [3lj . Only some nu- 
merical simulations are available in higher dimensions, 
and have provided the width distribution [U H|[ , which 
is a transient quantity therefore beyond the scope of this 
work. 

Over the last decade, the nonperturbative rcnormal- 
ization group (NPRG) has emerged as a powerful tool to 
investigate strongly correlated systems and critical phe- 
nomena (34)]. Besides a widespread application for equi- 
librium problems, the NPRG was further extended to 
study stationary (35l . |3(I and also transient nonequilib- 
rium systems (see e.g. ^33, In Ref. [H, the NPRG 

techniques were first applied to study the stationary KPZ 
equation in d dimensions. A simple ansatz was devised, 
referred to as leading-order (LO) in the following. This 
simple approximation already entails generic scaling and 
correctly captures the strong-coupling behavior in all di- 
mensions, providing critical exponents at the strong cou- 
pling fixed point [39( which are consistent with existing 
results in physical dimensions. 

A general and systematic framework to apply NPRG 
techniques to the KPZ equation is presented in Ref. (ioj . 
In this work, the symmetries of the KPZ equation are 
analyzed in detail and the related Ward identities are de- 
rived. A 'covariantization' associated with the Galilean 
invariance is introduced, which provides a simple geomet- 
rical interpretation of this symmetry, useful to formulate 
ansaetze respecting it. An appropriate NPRG approxi- 
mation scheme, strongly constrained by all these symme- 
tries, is devised. The second-order (SO) approximation 
(corresponding to a truncation at quadratic order in the 
response field), is explicitly implemented and the SO flow 
equations are solved in d = 1 to obtain the scaling func- 
tions and the amplitude ratio for the KPZ equation in 
this dimension [401 ] . The scaling functions turn out to 
compare very accurately with the exact ones Q, repro- 



ducing in particular all the fine structure of the tail on the 
correct scale. However, the one-dimensional case remains 
special due to the incidental existence of the generalized 
FDT and the NPRG flow equations in d = 1 are simpler. 
The numerical solution of the SO ansatz in arbitrary di- 
mensions is more involved. 

In the present work, we propose a simplification of the 
SO ansatz which allows us to solve the corresponding 
flow equations in arbitrary dimensions without cumber- 
some numerics. We will refer to this approximation as 
next-to- leading order (NLO). The simplification consists 
in neglecting the frequency dependence of some flowing 
functions in the integrals of the NPRG equations (see 
Sec. IIIip . This approximation drastically simplifies the 
frequency sector of the flow equations, enabling one to 
calculate analytically all the frequency integrals. It re- 
sults in a reduction of the computing time for the numeri- 
cal integration of the flow equations by a factor of around 
200, which in turn allows us to achieve an extensive study 
of the NLO flow equations in various dimensions. 

The physical quantities of interest in this paper are the 
correlation and response functions of the stationary KPZ 
equation in the strong coupling regime. From the correla- 
tion and response functions we extract the related scaling 
functions and determine universal amplitude ratios. We 
first study the one-dimensional case. The scaling func- 
tions obtained at NLO, although less precise than with 
the complete SO approximation, still compare accurately 
with their exact counterparts. We then investigate higher 
dimensions. We provide the scaling functions and univer- 
sal amplitude ratios in d = 2 and 3. To the best of our 
knowledge, no other comparable results are available in 
the literature. The obtained scaling functions show an 
exponential decay, but on a different scale than that pre- 
dicted by SCE and MCT QJ, |3, as in d = 1. We also 
determine the deviations in higher dimensions from the 
generalized FDT, which is only exactly fulfilled in d = 1. 

The validity of the NLO approximation can be assessed 
in two respects: in d = 1 by confronting the scaling 
functions with the exact results as mentioned before; in 
d > 1 through a detailed analysis of the critical expo- 
nents. First we find at NLO values in d = 2 and 3 that 
are consistent with results from numerical works reported 
in the literature [25|,[4l|. Moreover, we investigate the de- 
pendence of our results on the NPRG regulator, which 
provides us with some further check of the quality of 
the approximation. Although the NLO approximation 
deteriorates with the dimension, we argue that our NLO 
ansatz is reliable up to d ~ 3.5. The NLO ansatz is there- 
fore applicable to study the KPZ equation in dimension 
2 and 3, but cannot contribute to probe the existence of 
an upper critical dimension. 

The remainder of the paper is organized as follows. In 
Section [XXI we briefly review the nonperturbative renor- 
malization group formalism for out-of-equilibrium prob- 
lems. Section IlIII presents the symmetries and a general 
NPRG approximation scheme for the KPZ problem. We 
also specify the NLO approximation and derive the cor- 
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responding flow equations. The calculation of the fixed 
point properties is presented in Section [IV] We provide 
critical exponents, amplitude ratios and universal scaling 
functions. SectionfVl contains a summary and concluding 
remarks. Technical details are postponed in the Appen- 
dices. 



II. THE NONPERTURBATIVE 
RENORMALIZATION GROUP 

The general NPRG formalism for noneq uilibrium sys- 
tems is expounded in details in Ref . [36l . l3q ] . Its applica- 
tion to the KPZ equation is further specified in Ref. (ioj . 
Consequently, we recall for completeness only the main 
elements of this formalism, essentially following Ref. [4(| . 
In particular, the Fourier conventions used throughout 
this work are: 

f(uj,p)=[d d xdtf(t,x)e- t P- s+lult = [f(x)e- i ^ +iu ' t , 



(5a) 



d d p duj 
(2n) d 2tt 



ip-x — iu)t 



= f(p)e 



ip-x—iujt 



(5b) 



where x = (t, x) and p = (lo,p) and the factor (2-7r) d+1 is 
implicitly included in the momentum integral. 

The starting point is the field theory associated with 
Eq. ([!} , which can be derived from the standard Janssen- 
de Dominicis procedure {42| relying on the introduction 
of a Martin-Siggia-Rose [43| response field h and sources 
(j, j) . The generating functional reads: 

Z\j,j]=jv[h,ih] exp(-S[h,h) + J (jh + jh)), (6a) 

S[h,h]=J j/i (d t h -v\J 2 h- ~ (V/i) 2 ^ - Dh 2 ^ . (6b) 

Averages of an observable O over the disorder are ex- 
pressed in terms of the functional integral 



(0) = 



JV[h,ih}Oe- s V 1 '^ 
jV[h,ih]e- s \- h M 



(7) 



The NPRG formalism consists in building a sequence 
of scale-dependent effective models such that fluctuations 
are smoothly averaged as the (momentum) scale k is low- 
ered from the microscopic scale A, where no fluctuations 
are yet included, to the macroscopic scale k = 0, where 
all fluctuations are summed over [34|, [44|. For out-of- 
equilibrium problems, one formally proceeds as in equi- 
librium, but with the presence of the additional response 
fields and with special care required to deal with the con- 
sequences of Ito's discretization and with causality issues, 
as stressed in detail in [36|, [45[ . 

To achieve the separation of fluctuation modes within 
the NPRG procedure, one adds to the original action S 



a momentum and scale dependent mass-like term: 



(8) 



where the indices i,j € {1,2} label the field and re- 
sponse field, respectively h\ = h, hi = h, and summation 
over repeated indices is implicit. The matrix elements 
[R K (q)]ij are proportional to a cutoff function r(q 2 /K 2 ), 
with q = \q\, which ensures the selection of fluctuation 
modes: r(x) is required to vanish as x > 1 such that 
the fluctuation modes hi(q > k) are unaffected by A6> K , 
and to be large when x < 1 such that the other modes 
^ K )) are essentially frozen. Furthermore, AiS K 
must preserve all the symmetries of the problem and 
causality properties. An appropriate choice is 



R K (uj,q)=R K (q): 



v K q z 
2 -2D, 



(9) 



where the running coefficients v K and D K7 defined later 
(Eq. (|36D ). arc introduced in the regulator for conve- 
nience [391. Here we choose the cutoff function to be 



r(x) — a/(exp(x) — 1), 



(10) 



where a is a free parameter. 

In presence of the mass term AS K , the generating func- 
tional (|6aj) becomes scale dependent 



ZS~3\=\ v l h , ih] exp ( -S - AS K + / (jh + jh 



(11) 

Field expectation values in the presence of the external 
sources j and j are obtained from the functional W K = 
log-Z K as 

¥>(x) = (Mx)) = |^ , 0(x) = (/>(x)) = 5^r. (12) 



The effective action T K [ip, (p\ is now given by the Legendre 
transform of W K (up to a term proportional to R K ) [HI, 
11: 



T K [ip, <p) + W K \j,j] = J ji<Pi ~\ J Vi [Rk]ij ifij. (13) 
From T n , one can derive 2-point functions 



I r l 2) ] l li2(xi,X2,^,(p) 



6 2 T K [ip,<p] 
6ip il (x.i)8<pi 3 (x.2)' 



(14) 



and more generally n-point 1-PI (one particle irreducible) 
functions, that we write here in a 2 x 2 matrix form (omit- 
ting the dependence on the running scale K when n > 2) 



s n - 2 r {2 \ Xu x 2 ,<p^) 



5ip l3 (x. 3 )...6ip in (x„) 
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The matrix conventions are taken from Rcf. |36|. The 
exact flow for T K [ip, tp] is given by Wetterich's equation, 
which is in Fourier space |34[ 



d K T K = I Tr 



where 



G K = 



,(2) 



(16) 



(17) 



is the full renormalized propagator of the theory. When 
k flows from A to 0, T K interpolates between the micro- 
scopic model T k= a = S and the full effective action T k= q 
that encompasses all the macroscopic properties of the 
system [3(|. Differentiating Eq. ([To]) twice with respect 
to the fields and evaluating it in a uniform and station- 
ary field configuration (since the model is analyzed in its 
long time and large distance regime where it is transla- 
tionally invariant in space and time) one obtains the flow 
equation for the 2-point functions 

d K [ I?) ] 4 , (p) = Tr^ d K R (q) • G K (q) ■ (- \ T<§ (p, p, q) 

+lf >(p, q) • G K (p + q) • if >(-p, p + q)) • G„(q),(18) 

where the background field dependences have been 
omitted, as well as the last argument of the T^ n \ which 
is determined by frequency and momentum conservation 

Solving Eq. ([To]) (or Eq. (fI5|l) is in principle equivalent 
to solving the model. In practice this resolution cannot 
be performed exactly since these equations are nonlinear 
integral partial differential functional equations. Hence 
one has to devise an approximation scheme, adapted to 
the specific model under study and in particular to its 
symmetries. 



III. APPROXIMATION SCHEME 
A. Symmetries 

The approximation scheme must preserve the symme- 
tries of the problem. The KPZ action (|6b[) possesses three 
symmetries, besides the translational and rotational in- 
variances: i) the Galilean invariance, which translates 
for an interface as the invariance under an infinitesimal 
tilting, ii) the invariance under a constant shift of the h 
field, and iii) an additional time reversal symmetry only 
realized in d — 1. These symmetries can be expressed as 
an invariance of the action (|6bl) under the following field 
transformations 



h'(t, x) = x ■ v + h(t, x + Xvt) 
h'(t,x) = h(t,x + Xvt) 

ii) h'{t,x) = h(t,x) + c, 



(19) 
(20) 



where v and c are arbitrary constant quantities. As for 
the time reversa l sy mmetry, it can be encoded in the field 
transformation [la, [461 



h'(t,x) = -h(-t,x) 

h'(t, x) = h(-t, x) + %V 2 h{-t, x) 



(21) 



These symmetries are analyzed in details in [16|, |40| . 
where the corresponding Ward identities are derived. In 
fact, the symmetries i) and ii) admits a time-gauged 
form [47| , expressed as the following infinitesimal field 
transformations: 

. ;) f h'(t, x)=x- d t v(t) + h(t, x + Xv(t)) 



(23) 



h'(t,x) = h(t,x + \v{t)) 
ii') h'{t,x) = h(t,x) + c(t). 



As the variations of the KPZ action under these time- 
gauged transformations are linear in the fields, they 
also entail simple Ward identities, which extend with a 
stronger contents the ones deduced from the non-gauged 
forms 0. 

Regarding the shift-gauged symmetry, we simply recall 
a useful property of the n-point vertex functions, denoted 
from now on as 



r (I,m) ( 



Xl, 



i X/ +m ), 



(24) 



which stands for the p(™ _ ' +m ) vertex involving I (respec- 
tively m) legs - derivatives of T K with respect to tp (re- 
spectively (p) - with the I first frequencies and momenta 
referring to the ip fields and the m last to the ip fields. 
The shift-gauged symmetry entails that the n-point ver- 
tex functions in Fourier space satisfy the following prop- 
erty: 

r£"' n )(wi,i5i = 0,..., Pm+n _i) = wxS m xS n x (25) 

for any m > 1. This means that, apart from the con- 
tribution of f <pdt<p to r^ 1 ' 1 ^, the vertices vanish upon 
setting the momentum of one of the tp to zero. 

As for the time-reversal symmetry, we just give the 
Ward identity which relates the 2-point vertex functions 
in d = 1 



2Kr^)( P ) = --% 2 r(°. 2 )( P ) ) 



D 



(26) 



which will be useful in the following. 

Finally, the Ward identities ensuing from the gauged 
Galilean symmetry, relating the momentum and fre- 
quency sectors of vertex functions with different number 
of legs bare complicated forms and can be found in Ref. 
[40j . Let us stress instead a simple geo metrical interpre- 
tation of this symmetry proposed in [40J . A function /(x) 
is defined as a scalar under the Galilean transformation 
(|22p if its infinitesimal variation under this transforma- 
tion is Sf(x) = Xv(t) ■ V/(x). This then implies that 
J d d xf{x) is invariant under a Galilean transformation. 

With this definition h and VjV,-/i are scalars, but dth is 
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not. However, a scalar can be built from it if one intro- 
duces a covariant form of the time derivative D t h, denned 



AMx) = d t Mx)--(VMx)) 2 . (27) 

These three basic scalars can be combined together 
through sums, products and spatial gradients V to con- 
struct an action manifestly invariant under the Galilean 
transformation. The time derivative of a scalar is not a 
scalar itself, but one can again build a scalar quantity by 
defining a general covariant time derivative: 



b t = d t - AVh(x) • V, 



(28) 



which preserves the scalar property. This construction 
will ground the derivation of our ansatz. 



B. Rescaling 

Before we proceed, note that the bare parameters v 
and D of the original KPZ action ()6b[) can be absorbed 
upon rescaling the fields and time as 




t -> t/i 



(29) 



Followingly, the rescaled theory depends on a single ef- 
fective coupling constant 



A 



\ 2 D 



1/2 



9b 



and the rescaled KPZ action reads 
S[h,h] = f [h (d t h-V 2 h-^(Vh) 2 



(30) 



(31) 



which amounts to setting D and v to unity and changing 
A to ^J~gl in the original theory. Accordingly, the covari- 
ant time derivatives Eq. (|28[) in the rescaled theory read 



D t h(x) = dth(x) - ^(VMx)) 2 , (32a) 



D t = d t - VflfcVM*) ■ V. 



(32b) 



For future use, we give the relation between the original 
(l.h.s.) and the rescaled (r.h.s.) n-point vertex functions 
in Fourier space, which is 



r<[.™>(p 1 ,...,p, +ro _ 1 )| orl 



ginal 



1 Pi+m 1 ) | rescaled 

(33) 



with p' = (uj/v,p). Unless stated otherwise, we work 
from now on with the rescaled theory. 



C. NLO ansatz 

Our aim is to compute the physical correlation and 
response functions in the scaling regime. Consequently, 
our ansatz is designed to provide an accurate descrip- 
tion of the 2-point vertex functions, with arbitrary mo- 
mentum and frequency dependences, while approximat- 
ing higher order vertices. To do so, one can resort in 
general in the NPRG framework to the Blaizot-Mendez- 
Wschebor (BMW) approximation scheme, which is pre- 
cisely devised to compute momentum dependent n-point 
functions and has been proved to yield very accurate re- 
sults for the O(N) model [48|. However, the symmetries 
of the KPZ equation prevent from a direct application of 
the BMW scheme. Nevertheless, in Ref. 0, the BMW 
strategy was adapted to cope with the KPZ symmetries 
and an approximation scheme for the KPZ problem was 
built. An explicit ansatz for T K was proposed, which we 
call second order (SO). It is truncated at quadratic order 
in the response field h, while preserving the complete mo- 
mentum and frequency dependence of the 2-point func- 
tions. The SO ansatz reads 

r«M= j {<pfji(-i%,-v 2 )D t <p-$fZ( r &i,-v !l )<p 
- vV^-AVv 2 )^ + ^(-D 2 ,-V 2 )VV] } , (34) 

where the three running functions X £ {A, v 1 D}, are 
assumed to admit a formal expansion in a series of their 
arguments 



fZ(-Dt,-V J )= ^ a mn 



(-A 2 r(-v 2 



(35) 



that is, all the D t operators act on the left of the V ones. 
At the microscopic scale k = A, the three functions are 
constants f£ = 1, and the (rescaled) KPZ action (|3Tj) is 
recovered. 



1. Running anomalous dimensions 

We introduce two scale dependent parameters D K and 
v K defined as 



D K 



/°(0,0) , ^#(0,0) 



(36) 



to allow the fields and the dynamical scaling of time to 
renormalize during the flow. In the rescaled theory, both 
parameters are dimensionless and at the bare scale, D\ = 
v\ = 1. On the other hand, the function is fixed to 
unity at vanishing arguments /„(0, 0) = 1 for all k by 
the shift-gauged symmetry (see Eq . |3U|) . We define two 
anomalous dimensions associated with the two running 
coefficients as 

Vk = -Kd K \nD K , rf K = -Kd K lnv K . (37) 

Their connection with the physical critical exponents is 
specified in Sec. IIV1 
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2. Running vertex functions 

The calculation of the 2-point functions with the 
ansatz (j3~4"|) is straightforward and yields 



o, 

iuf* {U3,p) 



(38a) 
(38b) 
(38c) 



where p = \p\, and noting that the actual dependence 
of the three running functions /* is on ui 2 and p 2 . The 
shift-gauged Ward identity Eq. (j2"5")) further imposes the 
identity 



f*(Lj,p = 0) = l 



for all «. 



(39) 



The expressions of the 3- and 4-point functions are 
lengthy and reported in Ref. [40| . One can check that 
they fulfill the Ward identities related to the three sym- 
metries from Sec. IIII Al see Ref. [ioj . 

Note that d = 1 is a special case since the additional 
time-reversal symmetry leads to further constraints for 
the running functions. First the Ward identity (|26|) im- 
poses 



/«(w.P) = /« (w,p) = /«(w,p) 



and followingly 



V K =V K = Vk- 



(40) 



(41) 



Next, to consistently preserve this identity within the 
SO approximation, the entire f^(oj,p) function has to be 
fixed to one in this dimension (see Appendix C) 



ti(",P) 



(42) 



such that one is left in d = 1 with a single independent 
function f K and a single anomalous dimension r\ K [40| . 

In this paper, we work with the SO ansatz but im- 
plementing an additional approximation, which we re- 
fer to as next-to- leading order (NLO). The motivation 

is that, while all the T^'™^ functions and thus the flow 
equations greatly simplify in one dimension thanks to the 
additional time reversal symmetry, they remain compli- 
cated in generic dimensions. The NLO approximation is 
performed in the flow equations for the 2-point functions 
given by Eq. (|18[) . It consists in neglecting in the integrals 
on the right-hand side (r.h.s.) of the flow equation (JTSJ) 
the frequency dependence of the three running functions 
f£ , hence only keeping in the integrands the explicit fre- 
quency dependence of Ik' 1 ^. The NLO approximation 
thus consists in performing the replacement 



fZ(u,p)^fZ(p) 



(43) 



in all n-point vertex functions on the r.h.s. of the flow 
equation (|18|) . and in the propagator G K . The substitu- 
tion Eq. (|4"3")l in the 3- and 4-point vertex functions yields 

(2 1) 

that they all vanish except T„ ' ' which becomes 

T^iuj^^PuA) VfKPi +Pa\)- (44) 

One can check that the NLO approximation preserves 
the KPZ symmetries. 

The NLO approximation resembles the simple approx- 
imation achieved in Ref. [39} , to which we refer as leading 
order (LO). Indeed, the LO approximation is recovered 
by setting the entire function f£(u),p) = 1 (not only at 
p = 0) in all dimensions in the NLO flow equations. 
In dimension one, the two approximations LO and NLO 
coincide as f£ = 1 according to Eq. (|42p - unless one 
chooses to relax this constraint at NLO, as detailed in 
Appendix C. 



D. Flow equations 

We now derive the flow equations for the three running 
functions f*,XE {A, v, D}. For this, we consider the 
flow equations of the 2-point vertex functions, given by 
Eq. (fT5|). The propagator G K (uj,q) is defined in Eq. (fTT]). 
Using the expressions of the 2-point functions Eqs. (|38j) in 
the NLO approximation, i.e. with the replacement (|43p . 
and of the regulator R K given by Eq. ([9]), we get 



1 



2k K (q) Y K (w,q) 
Y:(u,q) 



where 



k K (q) = fZ(q) + D K r(q 2 /K 2 ), 
l K (q) = q 2 (fn(<l)+v«r(q 2 /K 2 )), 
Y K (uj iq )^iujf^{q) + l K { q ), 
P K (u,,q) = (u>ti(q)) 2 + (l K (q)) 2 , 



(45) 



(46a) 
(46b) 
(46c) 
(46d) 



and q = \q\. The derivative of the regulator matrix 



is 



dR(a)-( ° 1 2d ^(l)\ 



(47) 



where we have defined 



X K r{y) , y = q 2 /K 2 , X&{D,v}, (48a) 
-X K (^r(y)+2yd y r(y)). (48b) 



Kd K S*(y) 



The expressions for the propagator (|45|) . for the unique 

("2 r) 

non-vanishing vertex T K ' (|44|) and for the derivative 
of the regulator matrix (|4T|) are then substituted in the 
NPRG equation (fT8|) to get the flow equations for the 
(momentum and frequency dependent) 2-point functions 



and Tk'' z> (ui,p). The flow equations for the 



-.(0,2 



functions /* are deduced following Eq. 
obtain 



Wc finally 
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dj°{m,p) = 2g b f^p) 2 



(q 2 + (p-q)) 2 K(Q) 



d K fZ{w,p) = -2g/ K X (p) 



P K (u,q)d K S^(q)-2q 2 l K (q)k K (q)d K S:(q)j, 

-p- <iti(Q) UQ) P*{u, q) d,s°( q ) 

2P- qtiiQ) UQ) Uq) K{q) + (P 2 + P- q) tiiq) K(Q)(u> 2 fZ(q) 2 - U<?) 2 )] Q 2 d K S^(q) 

q 2 + (p- q) 



q 2 + (p- q) 



- Hp- qf^QY P K (u, q) d K S°(q) 



^P K (LJ,q^P K (n,Q) L 
to?- qtiiQf k K (q) + CO (p 2 +p- q) f^(q) 2 k K (Q)} q 2 l K (q) 8 K S v K (q) 



(49a) 



(49b) 



(49c) 



To shorten the notation, we introduced Q = \p + q\ and 
= w + uj. Notice that even if frequencies have been 
neglected in the / * functions within the integrands of the 
flow equations (|49| , these integrands still have an explicit 
(internal and external) frequency dependence, which is 
polynomial both in the numerators and denominators. 
This generates a non-trivial frequency dependence for 
the flowing functions /* on the l.h.s. of the NLO flow 
equations (|49|) . For consistency, the frequency indepen- 
dent flowing functions f*(p) on the r.h.s. of Eqs. (|4"9")l 
are obtained from the frequency dependent ones on the 
l.h.s. f*(za,p) evaluated at zero external frequency, that 
is f*(p) = f*(0,p). The integrals over the internal fre- 
quency w can be simply carried out analytically (see Ap- 
pendix A). 



E. Dimensionless flow equations 

As we intend to analyze fixed point properties, we in- 
troduce dimensionless and renormalized quantities (de- 
noted by a hat). Momentum and frequency are measured 
in units of the running cutoff K, 

p = p/n , w = w/(v k k), (50) 

and we define (dimensionless) renormalized functions as 

f$(w,$) = f${w,p)IX K , (51) 

using again the notation X to designate the three param- 
eters D, v 1 and A, with the conventions 

X e {D, v, A}, X K e{D K ,v K ,l}, ^6{CC0}- 

(52) 

D K and v K are the running coefficients Eq. (|36[) . which 
yields by definition 



We also introduce the dimensionless flow variable 

s = ln(/s/A) , d s = nd K , (54) 

often referred to as RG "time". The initial scale k = A 
corresponds to s — and the macroscopic scale k — > is 
obtained in the limit s — > — oo. 

The flow equations for the dimensionless functions are 
simply deduced from Eqs. (|50I51[) as 



dsf* = k8 k 



xJ K 



(55) 



with the notation (f52[) . This yields 

d.ft(w,p) = rfift(w,p) + (2 - tQw d*f*(m,p) 

+ pdpf*(zb,p)+iZ(zb,p). (56) 



The first three terms of Eq. ([56]) stems from the dimen- 
sional part. The nontrivial nonlinear contribution is cap- 
tured by the dynamical part encoded in the dimensionless 
integrals 



(57) 



The integrals I£ are given by the r.h.s. of Eqs. (|49j) 
where the substitutions for dimensionless quantities Eqs. 
(|5QI5ip are performed. Consequently, the bare coupling 
g b which appears linearly in the three integrals is changed 
to 



g b -> g b K d D K /v\ 



9 k, 



(58) 



which is dimensionless. In fact, when going to the di- 
mensionless renormalized quantities, all the dependence 
in the running coefficients D K and v K is absorbed into the 
running parameter g K , which is the only remaining inde- 
pendent coupling. Its flow equation is simply deduced 
from Eqs. ff3T|58|) as 



f*(m = 0,p=0) = l for all k. 



(53) 



d s g K =9 K (d-2 + 3rt v K -r,°), 



(59) 
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that is g K evolves only according to its dimensional flow. 

Following the definitions (|36I37[) , the running anoma- 
lous dimensions r\ v K and r\° are obtained from the zero 
momentum and frequency sector. Evaluating Eq. (156 
for X = v and X = Da,trb=p = 0, yields 



?7 °+/°(0 ! 0)=0, 
<+#(0,0) = 0, 



(60a) 
(60b) 



as /£(0,0) = f° (0,0) = 1 according to Eq. (f53j) . 

The two integrals in Eqs. (pOj) both have a linear 
dependence in the two anomalous dimensions according 
to Eqs. (|48I49I) . which we render explicit by writing 

t° (o, o) = ir c + iTnl + #° , (61a) 

£ (0, 0) = #X + « + (61b) 

The expressions of the J** terms are given for com- 
pleteness in Appendix A. Solving Eqs. (|60I61[) for the 
anomalous dimensions then yields: 



'7* = 



(i+/°°)(i+/-)-/r^ D 
(i+/°°)(i+jr)-/r^' 



(62a) 
(62b) 



The details to carry out the numerical integration of the 
flow equations (|56l59p are summarized in Appendix A. 



IV. FIXED-POINT PROPERTIES 



A. KPZ phase diagram 



We integrate the flow equations (|56I59[) in various di- 
mensions and for different initial values of the single di- 
mensionless bare coupling g\ = ^ d2 gb to determine the 
phase diagram of the KPZ equation, already presented 
in [lil, Let us summarize the result. In d < 2, the 
flow always tends to a nontrivial strong-coupling (SC) 
and fully attractive fixed point with gl c 7^ 0, and thus 
the interface always roughens. In dimensions d > 2, two 
different regimes can be reached, depending on the initial 
bare coupling g\. Below a critical initial value g A , the 
running coupling g K flows to zero and the Gaussian fixed 
point is reached. The corresponding interface is smooth 
and characterized by the Edwards- Wilkinson exponents 
z = 2 and \ = (2 — d)/2. Above g%, the running coupling 
g K flows to a strong-coupling fixed-point with gl° 3> 1. 
It describes a rough interface with x > and where the 
exponent relation Eq. (|68|) is fulfilled. The critical value 
3 A separates the basins of attraction of these two fixed 
points. Right at g\ = g%, the flow leads to an unstable 
fixed-point with < g^ T < gl° which drives the rough- 
ening transition (RT). 

In the following, we show how physical observables, 
which characterize the stationary KPZ growth, can be 



obtained from the fixed-point solution of the NPRG flow 
equations. We present our results for critical exponents, 
correlation and response functions and the associated 
universal scaling functions and amplitude ratios in di- 
mensions d = 1, 2 and 3. Since we consider only the 
strong-coupling behavior, we drop the explicit 'SC la- 
bel in the following, keeping in mind that all fixed point 
quantities, denoted by a star, are obtained at the strong 
coupling fixed point. 



B. Critical exponents 

In this section, we discuss the critical exponents ob- 
tained within the LO and NLO approximations in dimen- 
sions d = 1, 2 and 3. We first establish the connection 
between the anomalous dimensions r\° and r\ v K and the 
roughness and dynamical critical exponents x an d z - At 
a fixed point, the running anomalous dimensions attain 
their constant fixed point values r\° — > rj° and rf K 
From Eq. (|37p . we deduce that the running coefficients 
then acquire power law behaviors 



D K = D k 



(63) 



where Dq and i/q are two nonuniversal constants. On the 
other hand, the physical critical exponents x an d z are 
defined as the anomalous scaling of the frequency and of 
the correlation function as 



W~K 2 , C(t,x) - K~ 2X . 



(64) 



From Eqs. (|50I51[) , the dimensions of the frequency and 
of the correlation function in the NPRG framework are 

[lu] = v k k 2 , [C(t,x)] =K d - 2 D K /is K . (65) 

Comparing Eqs. (|63I64I65[) then yields 

z = 2-€ , X=(2-d + ^-r£)/2. (66) 

On the other hand, the flow equation (f5"9"|) of the coupling 
g K implies that any non-zero fixed point with g* =/= 
satisfies the relation 



(d-2 + 3^-r?°)=0. 



(67) 



Combining Eqs. (|67I66[) then enforces, as expected, the 
exponents identity 



z + X = 2 



(68) 



at any non-gaussian fixed point. 

We hence restrict the following discussion of critical 
exponents to the values of \- I n dimension d = 1, the 
two anomalous dimensions are equal and the exact value 
X = 1/2 is recovered. This result holds both at LO and at 
NLO, since both approximations coincide in d = 1 once 
the time reversal symmetry is satisfied. In dimension d = 
2, the NLO approximation yields x — 0.373(1), which is 
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« — >d= 1 NLO (FDT) 

d = 1 NLO (no FDT) 
o-^>d = 2LO 
o — °d = 1 NLO 
^-vrf=3LO 
d = 3 NLO 
°d=4LO 
<2 = 4 NLO 
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FIG. 1. (Color online) Variations of the roughness exponent \ 
with the cutoff parameter a from the NPRG approach within 
the LO and NLO approximations in integer dimensions d. In 
d = 1 the curves are computed with the NLO ansatz, either 
imposing the time reversal Ward identity Eq. (I26p (FDT) - in 
which case NLO coincides with LO - or relaxing it (no FDT). 



in better agreement with the numerical results than the 
X value obtained at LO (see table [TJ) . In dimension d = 
3 however, the value of \ is not significantly improved 
from LO to NLO and the discrepancy with the numerics 
remains noticeable. The results for the exponents are 
summarized in table HI 



TABLE I. Critical roughness exponent \ from the NPRG ap- 
proach at LO and NLO approximations in integer dimensions 
d. The values 'Lit.' report results from other works: the ex- 
act result in d = 1, and in higher dimensions, it corresponds 
to an average (with standard deviation) over values obtained 
in [H EH, IMQ. The error bars for the calculated NPRG 
values reflect the variations around the PMS values (see text). 



d 


1 


2 


3 


4 


X(LO) 


1/2 


0.330(8) 


0.173(5) 


0.075(4) 


X (NLO) 


1/2 


0.373(1) 


0.179(4) 




X (Lit.) 


1/2 


0.379(15) 


0.300(12) 


0.246(7) 



The variation of \ as a function of the cutoff param- 
eter a for different dimensions and approximations are 
displayed in Fig. [T] (note the discussion on the role of 
a in appendix A). The values of x m table U represent 
the minima of the x( a ) curves of Fig. [TJ The errors re- 
flect the deviations above these values, when a is varied 
within a total range of width 20 (d = 1), 10 (d = 2), 
1 (d = 3) and 0.5 (d = 4) around the minimum values 
apMS- To estimate the accuracy of the NLO approxi- 
mation, let us further analyze the x( a ) curves of Fig. [1] 
Whereas these curves appear essentially flat for dimen- 
sions one and two, their a dependence increases with 



the dimensionality. For large values of a we find that x 
roughly changes logarithmically with increasing slopes as 
d grows both at LO and NLO approximations (except for 
the d = 1 NLO with FDT curve, where x — 0-5 for all a). 
Moreover, for d = 2 and 3, the function x( a ) exhibits a 
stable minimum at both orders NLO and LO. For d > 3.5 
(curves in non- integer d not displayed) , a PMS keeps ex- 
isting at LO, whereas it disappears at NLO, which indi- 
cates a clear lack of convergence of our ansatz in higher 
dimensions. Consequently, the results are expected to be 
accurate for d — 1 and d = 2, seemingly less accurate but 
still reliable in d — 3, and not to be trusted for d > 3.5. 



C. Scaling form of the flowing functions 

We show in this section that the three functions f£ ac- 
quire a scaling form at the fixed point. In all dimensions, 
we observe that all the f£ functions take a nontrivial 
shape at the strong-coupling fixed point. They are bound 
to unity by definition (|53|) at vanishing momentum and 
frequency and their tails decay algebraically for large w 
and/or p. The fixed point functions f* are by definition 
solutions of the stationary equations d s f* = in Eqs. 
(|56|) . Besides, we verified that the decoupling property 
is satisfied. That is, we determined analytically the limit 
of the interaction integrals I£ (w,p) in the regime w 3> 1 
and/or p 3> 1 and checked that they all tend to zero. It 
follows that, at the fixed point, Eqs. (|56| reduce in the 
regime vo and/or p 3> 1 to the homogeneous equations 

Vff* + p dpf? (w,p) + (2 - r,:)w cW* (w,p) = 0, 

(69) 

with again the conventions Eq. (|52j) . One can easily show 
that the general solutions of Eq. (|69|) take the scaling 
form 

f*(vu,p)=p-^ ( x (m/f), (70) 

where z is defined by Eq. (|66p . The explicit form of 
the scaling functions £ x cannot be determined from the 
homogeneous Eqs. ([70)1 . However, they can be extracted 
from the numerical solution of the flow equations (|56|) 
by tabulating the values p 71 * /f {w,p) against vojp z . As 
shown below, the scaling of the physical correlation and 
response functions emerges from the form Eq. (|70|) of the 
fixed-point solution. 

Prior to this, let us determine the asymptotics of the 
scaling functions. We denote r = w/p z the scaling func- 
tion argument. As we find that the fixed point functions 
f*(^iP) are regular for all w and p values, constraints 
for the limits of the scaling functions £ x can be deduced. 
At vanishing argument r — > 0, C x (0) has to be finite for 
the limit w — > at fixed p to exist and followingly 

for p > vo xjz . 



fX(m,p)~( x (0)p-<>* 
At infinite r, C x must behave as 

C X (t)~C£t-^ for t ^ oo 



(71) 



(72) 
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for the limit p —> at fixed w to exist and with some 
constant C^. Followingly 



f?(va,p)~am-^/ z for &>j5* 



(73) 



Let us recapitulate the explicit leading behaviors of the 
scaling functions in the limit r — > oo : 



C X (r) ~ = 1, 



(74a) 
(74b) 
(74c) 



where the last constant in Eq. (|74cj) is fixed due to the 
shift gauged symmetry Eq. (f39"l) . 



D. Correlation and response functions 

The physical (dimensionful) correlation and response 
functions can be reconstructed from the (dimensionful) 
2-point vertex functions (see Appendix B) in the limit 
k — y 0, which corresponds to the fixed point. Let us first 
express the dimensionful fixed point functions in terms 
of the dimcnsionlcss ones. Following the definitions Eqs. 
(|51l) . they are given by 



f?(w,p) = xJ*(zu,p) = x K f; 



VJ p 



(75) 



The limit n — > at fixed vu and p is precisely equivalent 
to the regime w 3> 1 and/or p> 1 where the decoupling 
occurs, and where f* scales according to Eq. (f70|) . More- 
over, when the running scale k tends to zero, the behavior 
of D K and v K is controlled by the fixed point, where ac- 
cording to Eq. (|63|) they become power laws (with the 
nonuniversal constants v$ and D ). Hence, the physical 
dimensionful functions f* can be expressed in terms of 
the fixed point scaling functions ( x Eqs. ([70|l as 



f?(*>,p) = ^C B (W(^ z )) 

pv, 

f:(™, P ) = ^c fa/fop*)) , 



(76a) 
(76b) 
(76c) 



where the dimcnsionlcss argument has been expressed 
using Eq. (j50[) in terms of the dimensionful variables as 



r = m/p z = m/(v p z ). 



(77) 



Finally, according to Eqs. (|38I66I) . the 2-point vertex 
functions write 



M) (w, p) = WP z (ir( x (t) + C (r)) , (78a) 

(78b) 



V\ '{^,p) = --—^-( (r) 



P 



d-2+3 X 



On the other hand, the correlation and response func- 
tions are related to the 2-point vertex functions at the 
fixed point via (see Appendix B): 



1(0,2; 



(w,p) 



G{w,p) 



{vo,p)\ 2 ' 

(1,1)/ -A 



i(l,l) 



(m,p)\< 



(79a) 



(79b) 



From these relations we deduce that the physical correla- 
tion and response functions take in the stationary regime 
the scaling forms 



pd+2+X „ 



1 



(80a) 
(80b) 



where the two scaling functions F and H are defined as 



F(r) 
H{r) 



(tC A (t)) 2 +C(t) 2 ' 



(81a) 
(81b) 



The asymptotics of the scaling functions F and H in 
the limit r — > 00 can be simply deduced from the asymp- 
totics of the £ x functions in Eq. (|74l) as 

(82a) 
(82b) 
(82c) 



F(t) ~ C r~^+^ z , 



B(#(t)) ~ C 



v -_-(2+x/z) 



At vanishing r one obtains 

F(0) = C D (0)/(C(0))\ 
9f(F(0)) = 0, 

»(ff(o)) = l/no). 



(83a) 
(83b) 
(83c) 



Note that the expressions (|80|) were derived within the 
rescaled theory, whereas the KPZ equation is usually 
studied in its original version with three parameters A, D 
and v. According to Eqs. (|29l79p . the relation between 
the original (l.h.s.) and the rescaled (r.h.s.) correlation 
and response functions is 

C{za,p, D, u, A) = -^C(m/v,p,l,l,^/gb) (84a) 

G{zu,p, D, u, A) = -G{w/v,p, 1, 1, y/gfi. (84b) 
v 

The correlation and response functions in the stationary 
regime of the original theory are hence given by 



C(w,p) = 



DDo 



v i v -i p d+2+ X 



F 



vv Q p- 



G(tu,p) = — —H 
vvq p z 



vv§ p z 



(85a) 
(85b) 
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Equations (|85[) prove that the physical correlation and 
response functions endow a scaling form in the stationary 
regime, which evidences generic scaling. Let us empha- 
size that we did not assume the existence of scaling, it 
naturally arises from the presence of the fixed point and 
the form of the solution of the flow equations starting 
from any reasonable microscopic initial condition. The 
scaling functions F and H are hence universal with re- 
spect to the initial conditions up to the nonuniversal nor- 
malizations. Nevertheless, they still depend on the renor- 
malization scheme through the parameters and D and 
on the bare parameters v and D. This dependence can 
be removed via an appropriate normalization procedure, 
which will be described in the following. Prior to this, 
let us determine universal amplitude ratios, which are 
independent of the choice of normalization. 



E. Amplitude ratios 

The scaling form of the correlation function in real 
space is usually expressed as Eq. ([3]) with the scaling 
function F(y), behaving as 



F(y) 




(86) 



where Fo and F^ arc constants [22J. One can build a 
ratio of these amplitudes which is universal, which we 
present now. 

Let us first express the F Q and F^ amplitudes in term 
of the correlation function in Fourier space C(w,p) that 
we have calculated. For this, we consider the truncated 
correlation function 



AC{t,x) = C(t,x) - C(0,0), 



(87) 




which has the same asymptotic behaviors as C(t,x), 
namely 

AC(t,x) = 

and which has a well-defined Fourier transform 
*CM- r±f$ (e-<-- -!) <*.,», 

(89) 

where C(w,p) is given by Eq. (|85aj) . 

We first analyze the limit x — > 0. Noting that C(w : p) 
in Eq. (|85a[) only depends on we first change to the 
variables r = w / \vvap z ) and then to the variables u = 
p{vv t) 1/z in Eq. (JS9j) to obtain 



AC(i,0) = 2- 
where Td{x) is defined as 



d d u cos(tm 2 ) — 1 



(27r) d 



,2 X +d 



(90) 



F(r). (91) 



The expression Eq. (|91|) can be rewritten as 

T d = I(2 X /z)u d ( X ), (92) 
with the definition 

roc dr 



I(x) 



r X F(t). 



(93) 



The remaining integral u d (x) can be performed analyti- 
cally, yielding 

r 00 cos(u z ) - 1 

u d (x)=v d du ul+2x 

= _Vd [ r(-2 X A)oos(27r/z) X ^2/3, X <1 
z \n/2 X = 2/3. 

(94) 

The factor Vd is related to the angular integration volume: 
v- 1 = 2 d - 1 n d / 2 T(d/2). (95) 



Comparing the relations Eqs. (|88I90[) . we deduce the F a 
amplitude 



l^l(pu ) 2x/z u d ( X )I(2 x /z). (96) 



We now consider the limit t — > 0. From Eqs. (|85al89l) . 
we obtain 



AC(0,x)=2^I(0)K d (x, X ), 
w 



with 



K d (x,x) 



d d p e^- x - 1 



(27r) d p 2 x+ d 



(97) 



(98) 



With the substitution p ■ x — vcosi?, the x dependence 
can be factored out 



K d (x,x) = \x\ 2x w d ( X ), 



(99) 



and the remaining integral w d (x) can be performed ana- 
lytically 



Vd-l 



dv 



2tt J vW J0 
2(2x+d) 7r d/2 r ( x + d / 2 )' 



ddsm d - 2 ti(e lvcos t> - l) 



(100) 



Comparing the relations Eqs. (|88I97[) . we deduce the Fo 
amplitude 



DD 

Fa =2 -I(0)w d ( X ). 



(101) 



With the amplitudes (|96|) and (jlOip . one can build a 
universal ratio, that we define as 



R 



F a 



(102) 
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Indeed, inserting the expressions of the two amplitudes 
in this definition, all nonuniversal factors cancel out to 
give 



R = 



u d ( X )I(2 X /z) 



(2s,K(k(x)U(o)) 2 / 2 



(103) 



This amplitude ratio is computed from our NPRG corre- 
lation functions performing numerically the two integrals 
1(0) and I(2x/z) (see appendix A). 

Our results for R are summarized in table ITT] In d = 1, 
another amplitude ratio g(0) has been considered in the 
exact calculation of Rcf. (4J, which is related to R via 
g(0) = 2 1 / 3 i? (see appendix C). The exact result for q(0 ) 
is the Baik-Rain constant, which is 1.150439.... |54|. 
From our R, we obtain in d = 1 for g(0) the values 1.23(1) 
at NLO (this work) and 1.19(1) at SO 0o[. The NPRG 
treatment hence tends to overestimate the universal am- 
plitude ratio in this dimension. The result is less accurate 
at NLO than at SO, as expected since NLO is a lower or- 
der approximation. The accuracy at NLO is of the same 
order as that of the MC approximation, which predicts 
g(0) =1.1137 0,111. To our knowledge, no predictions 
for the universal amplitude ratio R in d ^ 1 are available 
in the literature. 



TABLE II. Amplitude ratio R from the NPRG approach at 
NLO approximation in integer dimensions d. Error bars cor- 
respond to deviations when varying the cutoff parameter a, 
using the same intervals as in Sec. IIV Bl The exact value in 
d = 1 is obtained from the Baik-Rain constant [IIH and Eq. 
(EH1). To our knowledge no (approximate) results for R in 
d — 2 and 3 are published. 



d 


1 


2 


3 


R NLO 


0.977(1) 


0.940(2) 


0.951(4) 


R exact 


0.9131 







F. Scaling functions in d = 1 

In one dimension, we compare the scaling functions ob- 
tained from our two NPRG approximations, namely the 
NLO (this work) and the SO ansaetze |4(j , with the exact 
scaling functions obtained in Ref. , essentially in order 
to assess the quality of the NLO approximation. In Ref. 
0, three scaling functions /, / and / are defined, related 
one to another by the following Fourier transformations 



m 



- cos(rfc 3 / 2 )/(r), 

7T 

dk 

— cos(fcy)/(fc), 



(104a) 
(104b) 



and normalized using a specific criterion (see Appendix 
C). We use lower-case letters /, / and / to refer to the 
correspondingly normalized functions. In order to com- 
pute the numerical Fourier transformations, our function 



/ is smoothened by a third order fit E 
can be found in appendix A and Rcf. 



fO) - details 




FIG. 2. (Color online) Comparison of the one-dimensional 
scaling function f(k); exact result [j, NPRG at NLO (re- 
specting FDT, this work) and NPRG at SO [!}. The inset 
shows the stretched exponential behavior of the tail with the 
superimposed oscillations, developing on the same scale k 3 ^ 2 . 
Note the vertical scale: this behavior sets in with amplitudes 
below typically 10~ 6 . 



« — o exact 
o— oNPRG SO 
NPRG NLO 




FIG. 3. (Color online) Comparison of the one-dimensional 
scaling function f(y); exact result Q, NPRG at NLO (re- 
specting FDT, this work) and NPRG at SO 0. 

The scaling functions f(k) and f(y) obtained within 
NPRG are compared with the exact results in Fig. [5] and 
121 respectively. One observes that the overall agreement 
with the exact result is manifestly excellent, both within 
the NLO and SO approximations. Examining the de- 
tails, two slight differences between the NPRG scaling 
functions can be noticed. First, compared to the SO ap- 
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TABLE III. Characteristic quantities relative to the one- 
dimensional scaling function /: position of the first zero 
feo, coordinates of the negative dip (kd,fd), coefficient of the 
stretched exponential 60 and pulsation of the oscillations do 
of; from the exact result 0], NPRG at NLO (respecting FDT, 
this work) and NPRG at SO [4(|. The error bars reflect the 
deviations when a is varied within an interval of width 20 
around opms- 



quantity 


exact 


NLO 


SO 


ko 


4.36236 


4.01(5) 


4.60(6) 


kd 


4.79079 


4.89(4) 


5.14(6) 


fd 


-0.00120 


-0.0147(6) 


-0.0018(6) 


ao 


1 

2 




0.28(5) 


b 


1 
2 


0.134(2) 


0.49(1) 



proximation and the exact result, the origin of / at NLO 
is slightly shifted downwards in Fig. [3jSecond, the tail 
of the function / in the NLO approximation is less accu- 
rately reproduced than at SO. Indeed, nontrivial features 
of this tail are highlighted in 0] : / decays to zero with 
a stretched exponential tail and with superimposed tiny 
oscillations around zero, only apparent on a logarithmic 
scale. A heuristic fit of this behavior for k > 15 is given 
by 1: 

f{k) ~ 10.9 fc" 9/4 sin (a k 3/2 - 1.937) e~ b ° kV * . (105) 

Within NPRG, at both NLO and SO approximations, 
the decay follows a stretched exponential on the correct 
scale A: 3 / 2 , but with a less accurate coefficient bo at NLO. 
(see table UTTj ). Moreover, contrarily to the SO approxi- 
mation, the NLO one essentially misses the oscillations 
below a typical absolute / value of around 10~ 6 (see the 
inset of Fig. [5] and [HI.) Additional numbers to further 
characterize the function f(k) are listed in table UTTJ po- 
sition of the first zero fcrj, coordinates of the negative dip 
(kd-,fd), oscillation pulsation ao and decay coefficient bo- 
To summarize, the scaling functions in d = 1 obtained 
from the NLO and the SO ansaetze are very similar, with 
the SO approximation providing in general a more ac- 
curate description, closer to the exact result. This is 
expected since the NLO ansatz is obtained from a SO 
ansatz by truncating the frequency sector. However, the 
loss of precision remains small at the NLO approxima- 
tion, given the drastic simplification of the flow equations 
it entails. We can hence confidently restrict to the NLO 
approximation in higher dimensions. 

G. Scaling functions in d = 2 and 3 

The scaling functions F and H are constructed from 
our data as defined by Eqs. (|81"j) . They are related to the 
physical correlation and response functions via Eqs. ()85l) . 
The absolute normalizations of the scaling functions and 
their argument contain nonunivcrsal factors. In order to 



define fully universal scaling functions, one must rescale 
their abscissa and ordinates with an overall factor, fixed 
by adjusting the normalizations in an arbitrary conven- 
tional way. We here choose a different criterion than that 
of Ref. [4j used in the previous section. We introduce 
three parameters Co, Go and To as follows 

C(w,p) = Co d + 2+ F N ( —^-j r ) ,(106a) 

G(zj,p) = Go—^Hn(^^to), (106b) 
wo v \vvoV J 

(the subscript W denoting normalized functions). Com- 
paring Eqs. (|85ll06j) yields 

F n (t) = ±F(t/to) , H n (t) = ±H(t/t q ). 

Oo <J0 

(107) 

We then fix the three constants Co, Go and To by choosing 
the following normalization criteria 

r°° dr = 

F N (0)=H N (0) = 1 , / —F n (t) = 1, (108) 

Jo 7T 

which yields 

C = F(0), G =H{0), To = ~^~t~t ■ (109) 

b (0) 

The scaling function Fn is displayed in Fig. [4] and the 
real and imaginary parts of H jy are displayed in Fig. [5] 
This constitutes the main result of this paper. First, 




1000 



m/p z 

FIG. 4. (Color online) Normalized scaling function Fn (rel- 
ative to the correlation function) from NPRG at NLO in di- 
mensions 1, 2 and 3. 

let us emphasize that, as causality is preserved by our 
approximation scheme and as the regulator term ensures 
that the functions remain analytic in the complex upper 
half plane, the real and imaginary parts of the response 
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d= l,Re 
o--od= 1 , Im 

6- ad=2, Re 
o- -o = 2, Im 

d=3, Re 

7- vrf= 3, Im 




1000 




1000 



FIG. 5. (Color online) Real and imaginary parts of the nor- 
malized scaling function Hn (relative to the response func- 
tion) from NPRG at NLO in dimensions 1, 2 and 3. 



FIG. 6. (Color online) Ratio of the scaling function Fn and 
of the real part of the scaling function Hn in dimensions 1, 
2 and 3, which illustrates the departure from the generalized 
FDT. 



function are expected to satisfy the Kramers-Kroning re- 
lations 



»(ff(r)) 



IT 



dr 



, 3(ff(r')) 



1 

—V 

7T 



(110a) 



(110b) 



where V denotes the Cauchy principal value. We verified 
numerically that these relations are perfectly fulfilled. 

Several quantities can be studied to further character- 
ize these functions. The first one is the universal ampli- 
tude ratio R, which we already discussed in Sec. IIVEI 
The results are reported in table [TTJ A second quantity to 
measure is the departure from the generalized FDT, i.e. 
the ratio Fn /^R(Hn) as a function of r = va/p z , which is 
depicted in Fig. [5] This ratio is constant in d = 1 where 
the generalized FDT holds, that is F n /R(H n )(t) = 1. 
In dimensions 2 and 3, the departure from one increases 
as t grows since the powers of the algebraic decay of the 
two functions Fn and $1(Hn) differ in d ^ 1 according 
to Eqs. (PI) . 

Lastly, one can examine the structure of the tail of the 
function F defined, in analogy with the one-dimensional 
case Eq. (fTM)) . as 



F(k) 



dr 



cos(Tk z )F N (T). 



(Ill) 



The asymptotic behavior of F is determined by the pole 
structure of Fn- To analyze it, we model our data for 
Fn with the fitting procedure described in Appendix A. 
We obtain an analogous pole structure as in d = 1 at 
NLO. Namely, the singularity lying closer to the real axis 
is a pure imaginary one Z\. This entails for F(k) an 
exponential decay on the scale k z of the form 



F(k) ~ exp(-&x k z ) k 



oo, 



(112) 



and no oscillations (see Appendix A). However, as the 
absence of oscillations appeared as an artifact of the NLO 
approximation in d = 1, we cannot settle whether the 
oscillations would persist in higher d in a higher-order 
approximation [56| . The coefficient of the decay is b\ ~ 
0.61(1) in d = 2 and b x ~ 0.19(1) in d = 3. 

This result can be confronted with the prediction from 
both the SCE and the MC approximations [H,[27|. For 
instance, the function h(k) is studied within the MC 
approximation [l8| . It is related to our function by 
n(k) oc F{k 1 / Z ), and it is predicted to decay exponen- 
tially with the asymptotic form 



n{k) ~ (Bk) 



(d-l)/(2a) -IB/cl 1 /* 



k 



oo 



(113) 



with some constant B. It would mean for the function 
F(k) = n(k z ) a decay on the scale k, and not on the scale 
k z as found in our analysis Eq. (|112[) . We find the same 
discrepancy as for the one-dimensional case, where our 
result agrees with the exact result. 



V. SUMMARY 

In this paper, we analyze the strong-coupling sta- 
tionary regime of the KPZ equation using NPRG tech- 
niques. We work with a simplified version of the full 
SO (quadratic in the response field) ansatz derived in 
Ref. [40] ■ The simplification, that we call NLO, consists 
in performing an additional approximation in the fre- 
quency sector of the flow equations, namely neglecting 
some frequency dependence within the integrals of the 
flow equations, while preserving a nontrivial frequency 
dependence of the flowing functions. This simplification 
leads to a drastic reduction of the computing time for 
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the numerical integration of the flow equations, and docs 
not spoil the global quality of the SO approximation, as 
evidenced in the one-dimensional case. 

Indeed, we compute the scaling functions of the one- 
dimensional growth and confront them with the results 
obtained with the full SO ansatz [4(| and with the exact 
results of Ref. We show that the NLO approximation 
essentially preserves the excellent agreement with the ex- 
act results that was found at SO, reproducing most of the 
fine structure of these functions, and with a restricted loss 
of precision. The main contribution of this work is the 
calculation of universal quantities in dimension 2 and 3. 
We compute the correlation and the response functions 
and prove that they take a scaling form in the stationary 
regime. We provide the associated scaling functions, and 
we also calculate universal amplitude ratios. They con- 
stitute our most important results. Let us stress that this 
is, to our knowledge, the first predictions for full scaling 
functions and universal amplitude ratios in dimensions 
different than one. These predictions could be compared 
with results from future large scale numerical simulations 
in d = 2 and d = 3. 

The critical exponents we obtain at NLO compare ac- 
curately with numerics in d = 2 and reasonably in d = 3. 
However, the NLO approximation deteriorates as the di- 
mensionality grows. Wc provide hints that it becomes 
unreliable above d > 3.5, such that the quantitative de- 
scription of the strong-coupling fixed point properties in 
dimensions higher than 3 requires a more sophisticated 
approximation. While high-dimensional correlation and 
response functions might not be of direct physical rele- 
vance, a proper description of these dimensions appears 
essential in order to probe the existence of an upper crit- 
ical dimension and to settle on the long-lasting debate 
on its existence (see e.g. [H?} and references therein). 
The present ansatz could be improved in several ways. 
A first option would be to consider the fully frequency- 
dependent SO ansatz (j3"4"]l . which was solved only in d = 1 
so far |40j . Another option would be to enhance the re- 
sponse field tp sector, which is truncated at quadratic 
order in the SO ansatz. It is not a priori obvious which 
of these two options would be the most relevant in order 
to improve the high-dimensional behavior of the approx- 
imation. This is left for further investigation. 
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APPENDIX A: NUMERICAL 
IMPLEMENTATION 

Solution of the flow equations 

In this appendix we present our implementation to nu- 
merically solve the coupled flow equations (|56l59p . First, 
for d > 2, the d-dimensional integrals over the dimen- 
sionlcss internal momentum q, which have the general 
form 

m = J -0^9(\q\)F(\p + q\), (Al) 

can be simplified to 2-dimcnsional integrals. Using hy- 
perspherical coordinates, Eq. (|A1[) can be rewritten as 

/>oo 

I(p) = ^ ± dqq d - l g{q) 
27T Jo 

x / d6 sin d - 2 9 F( y / p 2 + q 2 + 2pqcos9),(A2) 
Jo 

involving only one radial and one angular integral. 

The flowing functions /* (m, p) are discretized onpxro 
grids in frequency m and momentum p = \p\ space. The 
grids have a linear spacing Ap and Am and a maximal 
range of ;5 max and w ma x- The NLO ansatz allows one to 
perform the frequency integrals analytically (see below). 
The momentum integrals are calculated numerically us- 
ing Simpson's rule. For d = 1, there is a single integral 
over q and for d > 2 there are two integrals over q and 9. 
Due to the presence of the d s R K term, the integrands de- 
crease exponentially with q, which enables one to safely 
cut the momentum integral at a upper finite limit. For 
momenta p+q > p m ax the functions f£ are extended out- 
side the grid using power law extrapolations. This corre- 
sponds to the expected asymptotics of the flowing func- 
tions, at least close to the fixed point. For momentum co- 
ordinates which do not fall onto mesh points, is inter- 
polated using cubic splines. The derivative terms pdp and 
md^j are computed using 5-point differences. We studied 
separately the influence of resolution (Ap, A9 and Am) 
and of mesh size (j5 max and n7 max ) on the results. In the 
frequency sector, typical grids range up to m max = 30 
with a spacing Am =1/4. The precision of the double 
numerical integral over the momentum is of order 1CU 4 
for the typical resolution Ap =1/8 and A9 = 7r/20 and 
an upper integration bound of p max /2 = 10. We checked 
that the variation of all physical quantities by changing 
the discretization stays of the order 10 -2 smaller than 
variations coming from the dependence in the cutoff pa- 
rameter a. 

For the propagation in renormalization time s, we use 
explicit Euler time stepping with a typical time step 
As = —2.5 x 10~ 5 . Starting at s = from the bare 
action (/^ = fj( = f£ = 1), the three functions /* are 
smoothly deformed from their flat initial shapes to ac- 
quire their fixed point profiles, typically after \s\ > 10. 
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As a convergence criterion we calculate the exponents 
from Eq. (|66[) and check the deviation from the exponent 
identity Eq. (|68[) . We typically record the fixed point 
functions at s = —20, where this deviation is smaller 
than 10~ 8 . We also have to specify the bare initial cou- 
pling. Since the running coupling constant becomes very 
large in higher dimensions, it is convenient to consider 
the reduced coupling constant g™ d = g K Vd, which stays 
of order one. A bare initial coupling g™ d typically rang- 
ing between 0.5 and 5 is sufficient to reach the strong 



coupling fixed point. We have checked that the results 
for universal quantities are independent of the initial val- 
ues, as expected. 

Let us now comment on the analytical integration over 
the frequency. As the function P K (uj,q) is quadratic in 
frequency, the integrands in Eqs. (|49j) are rational frac- 
tions in Cj with denominators of order 6 and numera- 
tors of order at most 3. The integral over the internal 
frequency Cj can hence be carried out analytically. For 
instance 



dto 



1 



oo 2^p K (i,g)p K (o,g) 

30 dCj Cj 2 
oo 2* P K (u),q)P K (£l,Q) 

did 1 

-oo ^ (P K (0J,q)) 2 P K (n,Q) 



(A3a) 



(A3b) 



Finally, we give for completeness the explicit expres- 
sions of the integrals at zero external momentum and 
frequency involved in the determination of the running 



(A3c) 



anomalous dimensions Eq. ([62 



JOQ 



K 4 



dqq 



dq q 



9 k- 



g K - 



TTVd+2 



d+3 r(q 2 )k K {q) 
d+5 r(q 2 )(h(<i)) 2 

mmUq)) 4 ' 

°j A d+5 {dM?)) 
dq — h X 

Hqk K (q)) 2 -2k K {q)l K (q) 
d+2 r{q 2 ) 



o dq ti(q)(Uq)) 3 



2f x K {qWM - q 2 d^{f x MUq)/q 2 



(A4a) 
(A4b) 

(A4c) 
(A4d) 



dq- 



f+ 3 r(q 2 ) 



2 Jo m)(l K (q)) 3 
(qft(q)dik K (q) - 2k K (q)(f^(q) + 



(A4e) 



i+5 



9kI dq- 



-jf^(q)k K (q)+^v d +2X 



mm.m 3 

(qf X {q)dgk K (q)- 2k K {q)Cf^) + qdgft(q)) 
- 2f X (q)(dMq))/q + qd d (ti(q)Uq)/q 2 



(A4f) 
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where 



L(q) = q 2 (.t(q) + r(q 2 )). 



singularity is a complex one zq = &o + ibo (and sym- 
metries), which coordinates are robust in successive or- 
(A5a) dej-g Q f the fits. This entails for the Fourier transform 
(A5b) the exponential asymptotic decay F(k) ~ exp(izo) oc 
cos(aofc) exp(— 6ofc) and thus for the function F 



Data analysis 

The physical dimensionful single argument scaling 
functions F(vj/p z ) and H(zu/p z ) (associated with the 
correlation and response functions) are constructed from 
the large p and/or w sector of the dimensionlcss two- 
argument fixed-point functions f*(m,p). In this work, 
gridded values between 18 < p < 19 and 28 < w < 29 
where selected to achieve the data collapse. This collapse 
is very accurate in all the dimensions and for all values 
of a studied. This was illustrated in details for the one- 
dimensional case in Ref. (40j |. However, in order to per- 
form the numerical Fourier transformations Eq. (|104[) to 
obtain / and /, small numerical noise superimposing the 
raw scaling function / causes numerical artifacts. For 
this reason, we model our scaling functions in terms of 
a family of fitting functions, as explained in Here 
we devise a generalization of this Pade-type fit, which 
especially reproduces the correct asymptotic decay Eq. 
(|82|) . At third order, this fit writes for a general scaling 
functions F: 



FfitW = 



2 , 4 , 6 \ {d+2+ X )/2z 

CO + C 2 T Z + C 4 T 4 + C 6 T° 



1 + ClT 2 + C 3 T 4 + C 5 T 6 + C 7 T S 



with the additional constraint 

4 (d+2+x)/2z 

c = i- 

U 



(A6) 
(A7) 



Generalization of this fitting functions to higher orders 
is straightforward. The family of fitting functions (|A6[) , 
with the appropriate computed exponents (given in table 
HI and Co = 1 fixed by the normalization criteria (jl08p . 
turns out to perfectly adjust to our data for Fn and 
5R(i?jv) in dimensions 2 and 3. We checked that the third 
order fit is sufficient to get a satisfying convergence. In 
fact, the same family of fits, with a global exponent one 
(pure rational fractions), also perfectly models the func- 
tion Q(_ffjv( T ))/ T (which is consistent with the Kramers- 
Kroning relations Eq. (|110b[l ). 

Moreover, the analytical form (|A6|) allows one to calcu- 
late the singularities of this function in the complex plane 
and determine the asymptotic behavior of its Fourier 
transform Eoj ]. defined by 



F(k) 



dr 



cos(rfc)F(T). 



(A8) 



Indeed, the decay of the tail of F is controlled by the 
singularity of F lying the closest to the real axis. For 
instance, in the one-dimensional case at SO, the closest 



F(k) = F{k z ) ~ cxp(-b £; z ) cos(a fc z ) k 



oo. 



(A9) 



If the closest singularity to the real axis is a pure imag- 
inary one z\ = ib\ (and symmetries), there is no oscil- 
lation superimposed on the exponential decay. We em- 
phasize that the exponential decay for F , and hence the 
exponential decay on the scale k z for F is rooted in the 
C°° nature of the scaling function F. Only the existence 
of a non-analyticity (e.g. a divergence at the origin) in F 
could drive an alternative behavior, which is precluded 
in our NLO ansatz. 

Let us make a last comment regarding the numerical 
calculation of the integrals I(x) defined in Eq. (|9"3"|) . which 
appear in the expression of the amplitude ratio. The pre- 
cision can be improved by exploiting the known asymp- 
totics of F. Indeed, as the function F decays as a power 
law according to Eq. (|82p. the contribution of the tail can 
be determined analytically after some sufficiently large 
crossover argument t c , and followingly 



I(x) 



F(tc)t c 



-(2 X +d-zx)/z 



2\ + d — ZX 



dTT X F(t) 



(AlO) 



where only the bulk integral is computed numerically us- 
ing a simple trapezoidal rule. 



Cutoff dependence 

For the numerical solution of the NPRG flow equa- 
tions, a value for the cutoff parameter a has to be chosen. 
We emphasize that, although the NPRG flow equation 
([T5| is exact and physical quantities are obtained in the 
k — y limit, where the regulator R K= o vanishes, approxi- 
mations introduce a spurious residual dependence on R K . 
After fixing the functional form of the regulator, this de- 
pendence can still be assessed by varying the prefactor 
a of the regulator and observing the change in the phys- 
ical quantities. This procedure provides us with some 
accuracy estimate and allows us to check the reliability 
of the approximation used. According to the principle 
of minimal sensitivity (PMS), the different values quoted 
in this work for physical observables correspond to ex- 
tremum values with respect to changes in a [Hs| . We 
further expect that the overall a dependence is weak for 
the approximation to be accurate. 
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APPENDIX B: CORRELATION AND RESPONSE 
FUNCTIONS 

In this appendix, we establish the relation between cor- 
relation and response functions on the one hand and 2- 
point vertex functions on the other hand. Conventions 
arc taken from Sec. ITT1 and Rcf. jsfj. The 2-point corre- 
lation function matrix W^ 2 ' is defined as 



W ( 2 )( 



<5j(x)5j(x') 
<5 2 W 



<5j(x)<5j(x') 
8 2 W 



<5j( x ) 5 j( x ') Sj(-x)Sj(x>) / j—]—o 

W (2 <°>(x,x') W( 14 )(x,x') \ 
W( 1 ^ 1 )(x,x') W(°' 2 )(x,x') I 



(Bl) 



In a uniform and stationary field configuration, the mo- 
mentum and frequency dependence of its Fourier trans- 
form simplifies due to translational invariance in space 
and time as 



W (2) (q,q') = W (2) (q) (27r) d+1 <$( d+1 >(q + q' 
where the matrix V\A 2 '(q) writes 



W (2) (q) 



W (2,o) (q) W (i,D(q) 
W^-q) W(°- 2 '(q) 



(B2) 



(B3) 



The matrix W^ 2 ) is obtained in the NPRG approach 
as the limit when k — > 0, i.e. at the fixed point, of the 
renormalized propagator G K = wji 



(2) 



W (2) = lim G K 



lim 

K->0 



if } + K 



A2) 



(B4) 



(2) 

As the 2-point matrix T K in Fourier space is also diag- 
onal in momentum, i.e. has the same form as Eq. (|B2I) . 
it can be simply inverted, which yields 



W (2) (q) = 



1 



rl°' 2) ( q ) rr ,i; (q) 



(1,1), 



rl M, (q)| 2 V r * M) (-q) -r^ UJ (q) 



(2,0), 



(B5) 



In real space, the correlation function we consider is 
defined by Eq. §3§ and the response function is introduced 

as 



G(t, x) 



S(h(t,x)) 



Sj '(0,0) 



(h(t,x)h(0,0)) c . (B6) 



Taking their Fourier transforms yields Eqs. (|79al79bl) . 



APPENDIX C: KPZ EQUATION IN D = 1 

In this appendix, we discuss the special case of the one- 
dimensional KPZ equation, which exhibits the additional 
time-reversal symmetry. As stressed in Sec. lIII C| the SO 
ansatz reduces in d = 1 to only one independent running 
function f K and one anomalous dimension rj K according 
to Eqs. (li0TH2l . 



Definition of the scaling functions in d — 1 

First, let us give the definitions and normalizations of 
the scaling functions calculated exactly in Ref . 0] , which 
are denoted with lower case letters and which serves as 
our reference. The real space correlation function con- 
sidered in [H is 

c(t,x) = ([h(t,x) - h(0,0) -t{d t h{t,x))f) = -2C(t,x), 

(CI) 

which differs from the correlation function C defined in 
Eq. ([3]) by a factor (—2). The scaling function associated 
with it is defined and normalized as 



(t, (2AW) 1/3 y) 



q(y) = lim 
^ ' (XNH/2) 2/3 



(C2) 



where N = aD/v is a normalization factor with an arbi- 
trary dimcnsionlcss number a. In fact, the second deriva- 
tive of this function is mainly considered in Ref. [3| 



f(y) = \g'\y). 



(C3) 



as well as two other scaling functions related by the fol- 
lowing Fourier transformations 



f(k) = 2 dy cos(ky)f(y), 



(C4a) 



/(t) = 2 / dk cos(/er).f(fc^), (C4b) 
Jo 



which can be inverted according to Eqs. (11041) . 
From Eqs. (fC2|) and (|C3)l one deduces 



,1/3 



(2A 2 iVt 2 ) 
f(y) = lim ^ 

t— >CO 



2N -9 x c(t,x)\ x={2X2Nt2) i /3y , (C5) 
and using Eq. (|C4|) we get 



-fc 7 / 2 

f{k) = lim — x 

" 2Nt (2X 2 Nt 2 ) 2/3 



t—¥oo 



°° dr . |3/2 , rk 3 / 2 k \ 

— cos(rfc 3 / 2 ) c ——, , (C6) 



7T 



* ' (2X 2 Nt 2 ) 1/3 



where the correlation function on the r.h.s. of Eq. (|C6I) is 
now the Fourier transform of c(i, x). Comparing the pre- 
vious equation with the definition for /(t) in Eq. (| 104|) fi- 
nally yields the relation between the Fourier transformed 
correlation function c(w,p) and the scaling function f(r) 
as 



-k J / 2 I k 3 / 2 

f ^ = A™ 2 5 /3A4/3AT5/3 t 7/3 C I T 



* ' {2X 2 Nt 2 ) 1/3 ) ' 
(C7) 
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Normalization in d = 1 

We implicitly work in the long time and long distance 
limit where scaling occurs. Our scaling function F is 
related to the correlation function c(w,p) by 



c(a7,p) = (-2) — -4— F { - 



v 2 Vq p 7 / 2 



3/2 



(C8) 



according to Eq. (|85ap and with the additional factor 
(-2) coming from Eq. (fUT]) . Substituting Eq. (fC8| into 
Eq. (|C7|t . and recalling that v = D a due to Eqs. <|40IG3[) . 
we obtain 



a 



1/2 



F((2a5*) 1/2 ' 



(C9) 



The pure numerical factor a is then fixed according to 
some normalization criterion, which was chosen in [4| as 



/(0) 



1, 



which induces 



r°° dr = 
2 -F(r) = 2I(0). 
Jo 



7T 



Inserting a in Eq. (|C9[) then finally yields 

1/2 



/(r) 



.9* 



1(0) 



F(2(3,/(0)) 1/2 - 



(CIO) 



(Cll) 



(C12) 



that we use in the one-dimensional section HVFI to com- 
pare our NPRG results with the exact scaling functions. 



Amplitude ratio in d = 1 

In Ref. [H the universal amplitude ratio is defined as 

/>oo 

5(0) = 4/ dyyf(y). (C13) 



The exact result for g(0) is the Baik-Rain constant and its 
numerical value is 5(0) = 1.15039... 0, [H3] ■ This quantity 
can be expressed in term of / as [40| 



ff(o) = ^r( \ 



drr 2 /3/( r ), 



(C14) 



once the two integrals are carried out analytically. 

Let us work out the relation between g(0) and the am- 
plitude ratio R defined in this work. The integrals Eqs. 
(|94|100|) simplify for d = 1 and x = 1/2 to 

Ul (l/2) = -r(l/3)/(27r), (C15a) 
^l(l/2) = -1/2, (C15b) 

so that the amplitude ratio R in Eq. (|103p reduces to 

r(l/3) 1(2/3) 



R = 



7r(^7(0) 4 )V3 



(C16) 



On the other hand, combining Eqs. (IC12IC14[) yields 
Q(0) 1/3 r(l/3) 7(2/3) 

■ 9(0) - 2 7T (5,J(0)4)V3' (C1?) 

from which one deduces 

R = .g(0)/2 1 / 3 . (C18) 



Note that Eq. (|C18[) can be used to calculate the am- 
plitude ratio g(0) in d = 1 even if ^ / 1/2, as R is 
well-defined is this case. This is relevant in the next sec- 
tion. 



Relaxed FDT constraint in d = 1 

The results obtained at NLO in one dimension are very 
satisfying and this gives confidence in this approximation. 
However, the one-dimensional case is special since the ex- 
istence of the additional time-reversal symmetry greatly 
simplifies the ansatz and the flow equations. Indeed, this 
symmetry imposes the identities = f° and = 1. 
Let us clarify the origin of this last relation. In generic 
dimensions, a generic function f x (-dj,p) is allowed. In 
d = 1, this function is related by the time-reversal sym- 
metry to some vertices that are neglected in our ansatz. 
As a consequence, our ansatz is only compatible with the 
time-reversal symmetry and FDT relations if one imposes 
/ = 1. Accordingly, one may wonder about the impor- 
tance of this constraint and what happens if one relaxes 
it. 

First, let us illustrate the role of this constraint. Once 
the identity f* = 1 is imposed, the other Ward identities 
associated with the time-reversal symmetry are automat- 
ically satisfied. Indeed, setting f£ — 1 in the integrands 
of the dimcnsionlcss flow equations (|56|) . the two inte- 
grals I°(m,p) and I^(w,p) in Eq. (fSTf become equal. 
Hence, starting from time-reversal symmetric initial con- 
ditions f£ = /jy, both flow equations remain identical 
and the time-reversal symmetry is preserved along the 
entire flow. On the other hand, if the identity f£ = l 
is not imposed, the dimensionless function f£ starts to 
flow, which breaks the equality between I V K and 1° and 
the two functions and f° become different. The FDT 
relation is then violated. 

In order to check the influence of the FDT constraint, 
we compare results obtained from the NLO ansatz in 
d = 1 with and without imposing it - that is either with 
= 1 fixed or with f£(w,p) flowing respectively. We 
ran both situations and now discuss the results we ob- 
tained. We first observe that the roughness exponent is 
no longer bound to the exact result \ = 1/2 at NLO 
without FDT, but the deviation remains small, we find 
X = 0.520(5). Note that the exponent identity Eq. ([68]) 
is still fulfilled as Galilean invariance is preserved. The 
overall cutoff dependence when varying a is weak as 
shown in Fig. [T] (no FDT). From Eq. (|102|) . which is valid 
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1 100 



m/ P z 

FIG. 7. (Color online) Comparison of one-dimensional scaling 
functions Fn obtained from different NPRG approximations: 
SO, NLO with FDT and NLO without FDT. The scaling 
function from the SO ansatz [io[ also respects FDT. Inset: 
violation of generalized FDT for the NLO ansatz without the 
explicit constraint f^ = l. 



in any dimension and for any reasonable values of x an d 
z, we can compute the universal amplitude ratio R. We 
find R = 0.871(2) (no FDT) compared to R = 0.977(1) 
(FDT) and R = 0.9131 exact result. We finally compare 
in Fig. [7] the explicit shapes of the scaling function Fn 
with and without FDT, which appear reasonably simi- 
lar. The inset illustrates the violation of FDT, which 
increases for large r as the asymptotics of the ratio is 
given by F n /U(H n )(t) ~ r -(3+x)/(4-x) an d x ^ i/ 2 
without FDT. 
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